Angular momentum of a strongly focussed Gaussian beam 
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A circularly polarized rotationally symmetric paraxial laser beams carries ±ft angular momentum 
per photon as spin. Focussing the beam with a rotationally symmetric lens cannot change this 
angular momentum flux, yet the focussed beam must have spin \Sz \ < fi per photon. The remainder 
of the original spin is converted to orbital angular momentum, manifesting itself as a longitudinal 
optical vortex at the focus. This demonstrates that optical orbital angular momentum can be 
generated by a rotationally symmetric optical system which preserves the total angular momentum 
of the beam. 

PACS numbers: 41.20. Jb,42.25.Bs,42.25. Ja 



INTRODUCTION 

That optical and other electromagnetic fields can carry 
angular momentum is a direct result of the fact that they 
can carry linear momentum. Since the linear momentum 
flux density of an electromagnetic field is E x H/c, giving 



— Re(E X H*) 



(1) 



for the time-averaged momentum flux density of a time- 
harmonic electromagnetic beam, it is tempting to write 



j = lRe(r X (E X H*)). 



(2) 



as the corresponding angular momentum flux density. 
However, since the conservation law for linear momen- 
tum only contains the time-derivative of the momentum 
density, and the divergence of the flux (or the volume 
integrals of these, noting that the volume integral of a 
divergence of a quantity can be written as the surface in- 
tegral of the quantity) , it is premature to identify this as 
the actual angular momentum flux density. If the above 
expression was in fact the correct angular momentum 
flux density, then the angular momentum of a circularly 
polarized plane would be zero. Since the correct classical 
angular momentum density must agree with the classical 
limit of the quantum angular momentum density, this 
must be incorrect. This was recognized long ago, and 
a separation of the angular momentum of the field can 
be separated into spin and orbital components yields a 
result in agreement with the quantum result Q, 

The total momentum or angular momentum flux 
through a plane can be unambiguously determined by 
integrating the above quantities over the plane; this is 
a straightforward application of the integral form of the 
conservation law for angular momentum, and is there- 
fore independent of the actual choice of expression for 
the angular momentum density. It can be noted that 
the actual momentum and angular momentum densities 
(as opposed to the flux densities given above) contain an 
extra 1/c factor, and can be integrated over the whole 



beam to give the total momentum and angular momen- 
tum content. However, for the case of an infinitely long 
beam, these quantities are infinite, and it is better to 
consider only the fluxes. 

Separating the total angular momentum flux of the 
beam into spin and orbital components gives 



J = L + S 



(3) 



where the total angular momentum flux J is the sum of 
an orbital term L and a spin term S. A similar rela- 
tionship exists for the angular momentum flux densities 
and densities at any point. The distinction is simply 
that the spin angular momentum density is independent 
of the choice of origin; that is, it is invariant with re- 
spect to translations of the coordinate system. Since the 
local conservation of angular momentum as the beam 
propagates in free space cannot depend on the choice of 
the origin of the coordinate system, both the orbital and 
spin components must be individually conserved, as well 
as the total angular momentum, since there would oth- 
erwise exist a choice-of-origin dependent torque exerted 
by the beam on free space. These quantities can only 
change if the fleld interacts with matter . 

In general, the separation into spin and orbital com- 
ponents is not quite as straightforward as one might 
hope One can, however, write expressions for 

the spin and orbital components. The Cartesian compo- 
nents of the time-averaged spin angular momentum flux 
density s are Q, 111 



tQlva{EyEl)/oj 
eolm{E^E*)/LU 



(4) 



where Ex,y,z are the Cartesian components of E, the 
complex vector amplitude. This result can also be 
written in terms of the Levi-Civita symbol as si = 
ieo£ijkEjE^/{2uj), where the expression is summed over 
repeated indices and the real part is taken. The orbital 
components are 0,13,13 



ieoEjir x y)E*/{2uj) 



(5) 
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These give simple results for paraxial beams. For exam- 
ple, a circularly polarized paraxial beam of power P and 
frequency lo has a spin angular momentum flux of ±P/w, 
depending on the handedness of the polarization 0| . This 
result is equivalent to the quantum mechanical result of 
±h per photon. If the beam has a uniform phase over 
a plane, or, more generally, phase that is rotationally 
symmetric about the beam axis, then the orbital angular 
momentum flux density is zero. 

We will consider such rotationally symmetric beams 
here, and present specific results for a Gaussian beam. 



ANGULAR MOMENTUM OF A FINITE BEAM 

The case of spin angular momentum flux equal to P/lu 
results only in the paraxial approximation, as it depends 
on Ez being zero. If we consider a beam of finite width 
in its focal plane, then the beam will spread through 
diffraction, and will, at a sufficiently large distance, be 
propagating in a purely radial direction. That is, for 
large r, we must have Er ~ 0. In this case, the electric 
field is purely tangential, and the spin angular momen- 
tum density in polar spherical coordinates is 



Sr = eolm{EgEl)/uj, 



(6) 



with the other vector components being zero. For a ro- 
tationally symmetric beam of the type we consider here, 
Sr will be independent of the azimuthal angle 0. 

Therefore, the maximum possible contribution to the 
total spin angular momentum, of which, by symmetry, 
only the z component is non-zero, is s^-cos^, where 9 is 
the angle measured from the z axis. Integrating this over 
the beam must result in 15*2 1 < h per photon. A conse- 
quence of this is that a focussed beam cannot be purely 
circularly polarized everywhere in a plane, including the 
focal plane. If the beam is rotationally symmetric and lo- 
cally circularly polarized in the far field, the beam cannot 
be partially plane polarized in the focal plane, either — 
therefore, the electric field must have an axial component 
{Ez ^0). A similar argument using linear momentum 
shows that plane polarized beams must have Ez =/= in 
the focal plane as well. 

As an example, we will find the spin carried by a fo- 
cussed Gaussian beam. A paraxial Gaussian beam has a 
(scalar) amplitude in the far field of 



U = Uq exp 



(7) 



where k is the wavenumber and wq is the paraxial beam 
waist radius j^, or in terms of the beam convergence 
angle given by tan^o = 2/(fcwo), 



For maximum possible spin, we have Sr = eoU'^/uj, and 
the total spin angular momentum of the beam, in units 
of h per photon, can be found by integrating over a hemi- 
sphere: 



Sz = A/P 



(9) 



where 



A= exp(-2tan2 6'/tan2 6'o)sin6'cos6'd6' (10) 

^0 



and 



P^ eiij>{-2ta-a^e/taii^9o)sm9de. (11) 

^0 

This can be readily integrated numerically, and the result 
for all practically realizable Gaussian beams is shown in 
figure n The qualitative behaviour seen here is expected 
since Sz — O.SP/w for a dipole radiation field 




Beam convergence angle e^^ (degrees) 



FIG. 1: Spin angular momentum in h per photon for a non- 
paraxial TEMoo Gaussian beam. The solid line shows the 
exact result, the dotted line shows the small-angle approxi- 
mation. 



For ( 
jiving 



<C 1, small angle approximations can be made. 



Sz^l- 9l/4. 



(12) 



U = Uo exp(- tan^ 9/ tan^ 6*0). 



(8) 



The relative error in the change in spin {1 — Sz) is less 
than 0.01 for beam convergence angles 0o < 6.25°, and 
less than 0.1 for 9o < 21°. 

Although the spin angular momentum flux of the beam 
is reduced by the beam being focussed, a lossless rotation- 
ally symmetric optical system cannot change the total 
angular momentum flux of an electromagnetic field j^. 
Therefore, there must be a corresponding increase in the 
orbital angular momentum. This is in remarkable con- 
trast to the usual methods of generating optical vortices 
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which employ astigmatic or cyhndrical lenses or holo- 
grams designed to break the rotational symmetry. The 
key difference is that orbital angular momentum genera- 
tion by focussing depends on the initial presence of spin 
angular momentum, whereas astigmatic systems do not. 

An interesting question that remains to be answered is 
in what way the focussed beam carries the orbital angu- 
lar momentum. This is best addressed by considering a 
rigorous electromagnetic model of the beam. 



Multipole expansion 

A time-harmonic electromagnetic beam can be repre- 
sented as a sum of of electric and magnetic multipoles: 



E(r) = X! X! ^nniMnm + 6„mN„ 



(13) 



where M„,„ and M„,„ are the TE and TM regular mul- 
tipole fields, or vector spherical wavefunctions |^. Not 
only are these wavefunctions a complete orthogonal set 
of divergence-free solutions of the vector Helmholtz equa- 
tion (and hence solutions to the Maxwell equations), they 
are also eigenfunctions of the angular momentum opera- 
tor J^, with eigenvalues [n(n-|- 1)]^/^, and Jz, with eigen- 
values m. The spin and orbital contributions to the an- 
glular momentum can be calculated from the expansion 
coefficients a„,„ and bnm |M ■ 

The only non-zero multipole coefficients for a left- 
circularly polarized rotationally symmetric beam are 
those with m — 1. Thus, the total angular momen- 
tum about the z axis is ti per photon. The multipole 
expansion coefficients for the beam can be determined 
by an overdetermined point-matching method . For a 
beam of finite width, it is found that the total spin is less 
than h per photon (spin calculated in this way exactly 
reproduces the curve in figure^. The remainder of the 
angular momentum is orbital. 

Since the multipole expansion of the beam is known, 
the fields can be calculated at any point in space. The 
components of the electric field in the focal plane are 
shown in figure |21 For a strongly focussed beam such as 
is shown in figureEl the longitudinal (ie z) component of 
the field is significant, with a magnitude of « 0.3 times 
the transverse components. All components of the elec- 
tric field show secondary diffraction rings (note that the 
radial dependence of multipole fields includes a spheri- 
cal Bessel function). The phases of the x and y compo- 
nents are uniform, except for an increment of tt between 
successive diffraction rings, and, due to the circular po- 
larization, differ by 7r/2 from each other. The phase of 
the z component, however, shows a clear azimuthal de- 
pendence identical to that seen 'm I — 1 paraxial vortex 
modes. Since this vortex behaviour is only possessed by 
the longitudinal component of the field, this can be called 
a longitudinal optical vortex. 



Calculation of the Poynting vector shows that there is 
indeed a transverse component, which, since its handed- 
ness is uniform, is responsible for the transport of the 
orbital angular momentum. 



(a) 




FIG. 2: Electric field components of a strongly focussed cir- 
cularly polarized Gaussian beam, with a convergence angle of 
45°. The X, y, and z components of the electric field in the 
focal plane are shown in (a), (b), and (c), while (d) shows 
phase contours (with a spacing of 27r/20) for the z compo- 
nent, showing azimuthal variation of phase as seen in vortex 
beams. 
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FIG. 3: Poynting vector of a strongly focussed circularly po- 
larized Gaussian beam, with a convergence angle of 45° . The 
transverse part (x and y components) of the Poynting vector 
in the focal plane are shown. 

As the beam is more strongly focussed, the magnitude 
of the longitudinal (z) component of the field increases, 
and the orbital angular momentum increases as a result. 
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The same increase can also be considered to result from 
the decrease of spin angular momentum, along with the 
conservation of total angular momentum. The change in 
the angular momentum and the growth of the longitu- 
dinal optical vortex is smooth and well-behaved as the 
convergence angle of the beam is increased, with no sud- 
den qualitative or quantitative changes. As the beam is 
more strongly focussed, the diffraction rings also become 
more prominent, but this does not affect the angular mo- 
mentum of the beam. 

DISCUSSION 
Torque density acting on lens 

In light of the above considerations, the transformation 
of spin angular momentum to orbital angular momentum 
by the lens must result in a reaction torque density about 
the z-axis acting on the lens that depends on the choice 
of origin. While such a torque acting on empty space is 
unacceptable on physical grounds, it is not only entirely 
reasonable, but expected, in the case of a lens. 

One need only to consider a simple ray picture of the 
action of a lens, not centered on the z axis, on a ray 
parallel to the z-axis. The focussed ray will generally 
not pass through the z-axis, and is not parallel to the z- 
axis, and hence carries orbital angular momentum about 
the z-axis. Consequently, there must be a reaction torque 
density acting on the lens. Note that this torque density 
depends on the choice of origin. The component about 
any axis parallel to the beam axis of the total torque 
acting on the lens is, of course, zero. 

Rotation in optical traps 

The presence of orbital angular momentum in the focal 
region suggests that orbital motion of absorbing or reflec- 
tive spherical particles should be observable in optical 
tweezers. However, since particles trapped in a Gaus- 
sian beam trap will be located on the beam axis, all that 
will be observed will be spinning of the particle about its 
axis. In order to observe this orbital angular momentum. 



it would be necessary to use a multi-ringed beam with 
zero orbital angular momentum, for example a Laguerre- 
Gauss mode LGpo, with p > 1, or a Bessel beam, as the 
input to the objective lens of the trap. In this case, par- 
ticles trapped in one of the rings rather than the central 
spot can be be expected to undergo orbital motion. 



CONCLUSION 

We note that focusing a circularly polarized beam pre- 
serves the total angular momentum flux of the beam 
about its axis. However, the spin component of the an- 
gular momentum flux is necessarily reduced as the beam 
is more strongly focussed. Due to the conservation of 
total angular momentum when the beam is focussed by 
a rotationally symmetric optical system, there must be 
a corresponding increase in the orbital angular momen- 
tum flux. This result is remarkable in that it predicts the 
generation of orbital angular momentum by a rotation- 
ally symmetric optical system, in apparent contradiction 
with common expectation. 

This orbital angular momentum is carried by the axial 
component of the electric field, which has the typical 
exp(i0) dependence of charge 1 optical vortices; we call 
this a longitudinal optical vortex. 
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